A two-step density-matrix renormalization-group study of coupled Luttinger liquids 
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We report a two-step density-matrix renormalization-group computation of the equal-time single- 
particle Green's function, the density-density correlations, and the low-frequency spectral weight 
function of a spinless fermion model in an anisotropic two-dimensional lattice at half-filling. We 
find that at weak couplings the density-density correlations have the universal decay of a Fermi 
liquid; the spectral weight function displays a sharp quasi-particle peak. But in the vicinity of a 
quantum critical point, these correlations strongly deviate from a Fermi liquid prediction and a 
pseudogap opens in the spectral weight function. 
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I. INTRODUCTION 

The metallic phase of interacting three-dimensional 
(3D) electron systems is described by the Fermi liquid 
theory (FLT) proposed by L. Landau 1 . The Fermi liquid 
(FL) is a phase of matter in which low energy charged 
excitations and the long-distance behavior of correlation 
functions are essentially similar to those of the non- 
interacting electron system. They carry spin and are 
described as weakly-interacting quasi-particles. A mi- 
croscopic justification of the FLT was given through the 
many-body perturbation theory^ 3 - and the renormaliza- 
tion group (RG) 4 ~ 6 . However neither the many-body per- 
turbation theory nor the RG provide a complete demon- 
stration of the emergence of a FL from an interacting 
electron model. The many-body perturbation theory ne- 
glects the competition between different channels and 
focuses only on the metallic phase. The RG analyses 
the flow of the interaction toward the FL fixed point, 
but it does not yield the quasi-particle spectra. Fur- 
thermore, both approaches are restricted to the weak- 
coupling regime. Indisputable FL behavior has been the- 
oretically shown only for impurity models^ or model in 
infinite dimensions 8 . Recent development into this diffi- 
cult problem involved the use of string theory 9 . 

Electron-electron interactions have a dramatic effect 
in the one-dimensional (ID) metallic phase. No matter 
how small, they completely destroy the quasi-particles. 
An alternative to the FLT for ID metals is Haldane'si^ 
Luttinger liquid theory (LLT). The central assumption of 
the LLT is that the low-energy excitations and the long- 
distance behavior of correlation functions of ID metals 
are similar to those of a model introduced by Luttinger 11 . 
These excitations are density fluctuations which propa- 
gate with different velocities for the spin and charge. In 
a Luttinger liquid (LL), unlike a FL, the decay of the 
correlation functions is non-universal. 

There is a significant interest in the question of the 
evolution of the Luttinger liquid when going from D = 1 
to D > 1. This is relevant to the physics of quasi-one- 
dimensional organic conductors 1 ^ where pressure or tem- 
paruture can induce a crossover from an LL to a FL or an 
ordered phase. The dimensional crossover has been stud- 
ied by various approaches. These include analytic con- 
tinuation from D=l to D = 1 + e 13 , perturbative renor- 



malization group (RG) on weakly-coupled LL 1 ^ or on 2D 
system with weak interaction 4 , functional integral^, and 
generalized dynamical mean-field theory (DMFT) 15 ' 16 . 
These studies conclude to a FL ground state in 2D. How- 
ever, Anderson and coworker s 1 i 18 have argued that, a 
different scenario due to strong interaction could take 
place. They pointed out that despite the RG being rel- 
evant, the resulting 2D system could nevertheless be a 
non-FL. The effects of the interactions could be so dra- 
matic that if the transverse hopping is not strong enough, 
the electrons would remain confined in the chains. Coher- 
ent quasi-particles would form only when the transverse 
hopping exceeds a treshold. This issue has recently been 
reexaminated in the framework of the functional RG 19 . 
A regime with confined coherence was predicted in the 
strong interaction regime. The non-FL mechanism sug- 
gested in Reb 17 ' 18 could occur for instance in the vicinity 
of a quantum critical point (QCP) where interaction ef- 
fects are very strong. 

In this paper, we use the two-step density-matrix 
renormalization group method 20 to study the possible 
emergence of FL and non-FL behaviors on an interacting 
electron model close to a QCP. Our results are consistent 
with a FL ground state in a 2D model for weak interac- 
tions. We also show that as a QCP is approached, the 
system enters a non-FL regime. This is captured by the 
behavior of the exponent K of the density-density corre- 
lation which shows a strong renormalization towards its 
FL value for V < 1 and which is only weakly renormal- 
ized in the vicinity of the ID quantum critical point. The 
evolution from a FL to a non-FL is also oberved in the 
low frequency spectral weight function. 



II. MODEL AND METHOD 

We concentrate on the following quasi-one-dimensional 
spinless fermion model on a finite lattice of size L x , L y 
in the x, y directions respectively: 

H = ~t x ^ { C t,iy C i* + hiy + k - C -) ~ M m.-iy ~ 

ty Y ( c L; y c w„+i + h.c.) + V Y n ix ,i y n ia+ i t i y . (1) 
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We are interested in the situation where the hopping pa- 
rameter t x along the x direction is far larger than the 
interchain hopping t y , t x 3> t y . The interaction V is cho- 
sen such that when t y = 0, we are in the LL phase, i.e., 
V < 2t x . We will restrict ourselves to case where the 
electron density is at half-filling, no — N e /(L x L y ) = 1/2, 
where N e is the total number of electrons. It has been 
shown in Ref. 20 that this type of anisotropic model 
may be studied using the density-matrix renormalization 
group (DMRG) method^. In this approach, the DMRG 
is applied in two steps. 

In the first step, we use the DMRG to construct an 
approximate, yet well controlled, low-energy Hamiltonian 
Ho,iy for an isolated chain i y Hamiltonian 

H 0liv = -t x J2( c La y c **+h*y +h.c.) - fiJ2 n **^y + 

(2) 

In order to allow interchain dynamics, Hqa is obtained 
by targeting the ground state of the nominal filling 
N ex /L x — 1/2, where N ex is the number of electrons 
on the chain i y . We also target ground states of N ex ± 1, 
N ex ± 2, N ex ± 3,... until the lowest state of a sector 
N ex ± k is higher than the highest state kept in the N ex 
sector. 

In the second step, the full 2D Hamiltonian (JTJ is pro- 
jected onto the basis constructed from the tensor product 
of the single-chain eigenfunctions; this projection yields 
an effective one-dimensional Hamiltonian for the 2D lat- 
tice, 

H « E Ha,i y - t y J2(4 v c ly+ i + h.c). (3) 

ly ly 

Ho,i y is diagonal, its element are the DMRG eigenvalues, 
cj , Ci y , and hi v are the renormalized operators in the 
single chain basis. These are vector operators made of 
local operators on each site of a chain i y . It is clear that 
during the passage from the first to the second step, this 
method is different from the conventional DMRG in that 
the truncation is not done through the reduced density 
matrix. The truncation is done rather like in the real 
space RG method^. But if t y remains small with respect 
to the energy width of the states kept, as in the Wilson 
approach for the Kondo problem^, this algorithm can 
retain high accuracy as we will show below. 

III. TEST ON THE NON-INTERACTING CASE 

Let us first analyze the performance of this DMRG 
algorithm for the case V = which enjoys an exact so- 
lution. The exact single particle energies and wave func- 
tions for open boundary conditions are respectively: 





16 x 17 


32 x 33 


64 x 65 


AE 


1.950841 


1.035854 


0.534903 


5e(t y = 0.05t x ) 


2.0 x 10~ 8 


< 10" 8 


2.4 x 10~ 7 


8e{t y = O.li*) 


1.0 x 10 -8 


4.6 x 10 -7 


2.9 x 10~ 5 



TABLE I: Energy width AE of the in 2 = 80 states kept and 
error 8e in the ground-energy per site as function of the lattice 
size. 



Ck; m = —2t x cosk Xl — 2t y cosk ym , (4) 



where k Xl — ln/(L x + 1), k ym = mn/(L y + 1), I = 
1, L x , m = 1, L y . The ground-state energy 

|k ! , m |<fc F 

the single particle Green's function between two points 
of coordinates (i x ,i y ) and (j Xl j y ), 

go{{ix,iy); (jx,jy)) =< cJ^CjWy >> ( 7 ) 

and the density-density correlation between these points 

Co((i x ,iy);(j x ,jy)) =< (n^ -no){n jx ,j v - no) >, (8) 
with n a7s , ay = cl^ ay c amtay , may be readily computed. 

9o((i X ,iy)\ Ux,jy)) = E ^kfe,%)-0k(ix,ij/), (9) 

k( ,m | <feF 

Co{{i x ,iy);(jx,jy)) is obtained from g ((i x , i y ); (j x , j y )) 
by the using Wick's theorem. 

In comparing the DMRG to this exact result, we em- 
phasize that although the exact solution is trivial in mo- 
mentum representation, for a real space method such as 
DMRG it remains a difficult challenge. However, unlike 
the exact solution, the DMRG can readily be extended 
to the case V ^ without difficulty. In the DMRG, we 
kept up to mi = 384 states during the first step with 
up to L x = 130. For this value of mi, the truncation 
error is virtually zero. Among the 4 x m\ states of the 
superblock, we kept a subset of up to nri2 = 108 states 
during the second step. These yield size dependent en- 
ergy widths AE — e,„ 2 — eo, where eo and e m2 are respec- 
tively the lowest state and the highest state kept in a 
chain i y (see Table [J). The key to retain accuracy during 
the second step is to choose t y such that, t y <C AE for 
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FIG. 1: Density correlation n(x) as function the position x for 
the middle chain of 2D systems: DMRG with V = 0, t y = 0.05 
(circles); the exact result with V = 0, t y = 0.05 (crosses); 
interacting 2D system with V = 1.5, t y = 0.1 (triangles); 
single chain with V = 1.5 (diamonds). 



a given L x . We show for instance in Table U the error in 
the ground state energies for t y = 0.05t x and t y = 0.1t x 
for L x x L y = 16 x 17, 32 x 33, 64 x 65. For t y = 0.05^, 
there is an excellent agreement with the exact energy for 
all sizes shown. Note that for the 16 x 17 and 32 x 33 
lattices, the limitation to only 8 digits is due to the fact 
that we set the error to 10~ 6 in the diagonalization of the 
Hamiltonian. We could easily reach a smaller error with- 
out significant additional work. The agreement remains 
excellent for t y = Q.\t x except for a 64 x 65 lattice. At 
this size, the difference between the DMRG and the exact 
energies is two orders of magnitude larger. The crossover 
temperature from ID to 2D is given by Tx ~ ty/ir—, for 
finite size systems in the ground state, this translates to 
6E w t y , where SE is the size-dependent lowest excitation 
in ID. We find that if SE/t y < 1 and AE/t y « 10, the 
accuracy is almost independent of the system size. For 
instance, when L x x L y = 32 x 33, 64 x 65, 130 x 131, 
we find respectively, SE = 0.0951, 0.0483, 0.0240. By 
respectively choosing t y — 0.1, 0.05, 0.025, we obtained 
6e < 10~ 6 . Hence, if we keep the ratios SE/t y and AE/t y 
constant, we can access the 2D regime in large systems 
while retaining very good accuracy. 

In order to limit the memory load, we computed the 
correlation functions in the central chain along the x di- 
rection, i y = L x /2 + 1 (longitudinal direction), and in 
the central chain along the y direction, i x — L x /2 (trans- 
verse direction). The local density no(x) (FigflJ, go(x) 
(Fig©, and C (x) (Fig© computed with the DMRG 
show a very good agreement with the exact result. We 
verified that the asymptotic behavior x 3> 1, go{x) ocl/x 
and Co (x) oc 1/x 2 of the exact result is satisfied by the 
DMRG. For instance in Figj2ja,b), for go(x) the largest 
difference Sgo between the DMRG and the exact result 
is seen in the tranverse direction at the largest distance 
x = 31 for which Sgo = 0.0002. The agreement is even 
better for Cq(x), SCo(x) rss 10~ 5 in the direction of the 
chains. For both the DMRG and the exact result, in 
the transverse direction, Cq(x) for x > 3 falls below 10~ 6 
which is the error set in the diagonalization of the Hamil- 
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FIG. 2: (a),(b), respectively go(x), Co(x) for a 64 x 65 sys- 
tem, for V — 0, t y — 0.05, the origin is on middle of the lat- 
tice, longitudinal direction: DMRG (circles), the exact result 
(crosses); transverse direction (only go(x) is shown): DMRG 
(squares), exact result (pluses). (c),(d) respectively g{x), 
C(x) for a 63 x 63 system, longitudinal direction: for V = 1.5, 
t y — (circles), for V = 1.5, t y = 0.1 (squares); transverse 
direction (only g(x) is shown): for V = 1.5, t y — 0.1 (dia- 
monds) . 

tonian. Hence, it was not shown. For this reason, we 
will exclusively concentrate on the correlation along the 
chains when analyzing the interacting case. 

IV. DENSITY-DENSITY CORRELATIONS 

For t y — and V ^ 0, there also exists an exact 
solution^. The model is in a LL phase for V < 2t x and in 
a charge density wave (CDW) phase for V > 2t x . In the 
LL phase, the asymptotic form of the Green's function is 
9id{x) ~ A g /x 1+a , where a is the anomalous exponent. 
The dominant two-particle correlations are the density- 
density, Cin(x) = A c /x 2K with a = {K+l/K)/2-l. At 
V = 2t x , there is a ID QCP. When open boundary con- 
ditions are applied, the sites at the ends generate strong 
Friedel oscillations. These oscillations decay very slowly 
from the ends and interfere with the normal density os- 
cillations Cid(x). This behavior of the ID system can be 
reproduced by the DMRG with extremely high accuracy. 

When t y sets in, it is expected that either the system 
will be dominated by the single particle correlation, hence 
the ground state is a FL, or the density correlations Cm 
would freeze yielding an ordered two-dimensional CDW 
state. We did not find any evidence of CDW long-range 
order (LRO) when we start from the disordered ID chain. 
It is important to note that the same method was used to 
study coupled spin chains and found LRQ^ii as expected. 
Let us further discuss the reliability of this result. Since 
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FIG. 3: C(x) for V = 1.5, t y = 0.054: (a) 66 x 67 (circles), 
102 x 103 (squares), and 130 x 131 (diamonds) systems; (b) 
L — 102 (ID, circles) and 102 x 103 (squares); the vertical 
dotted lines show the limit of the retained data used for the 
extrapolations. 



the DMRG is highly efficient in the interacting ID case, 
the level of accuracy in ID between the cases V = and 

V is comparable. Hence, studying the full Hamil- 
tonian [TJ which is no longer exactly solvable, with the 
DMRG is not more difficult than the V = case. The 
DMRG results of the 2D interacting case will be as good 
as those of the V = case, provided that AE does not 
decay sharply. When V ^ 0, an odd value of L x is cho- 
sen in order to not to frustrate the CDW correlations 
during the lattice growtb~. This choice also has the ad- 
vantage of showing a sharp contrast for the behavior of 
the Friedel oscillations between ID and 2D. We find for 
instance that in the LL phase for V — 1.5, L x = 63, 
AE increases from its V = value AE = 0.53 to 0.95. 
This implies that for the same value of t y , the accuracy 
would be better for the V = 1.5 case than for the V = 
case. However, for L x = 63 and V — 1.5, SE increases to 
0.11. For this reason, we have to increase the transverse 
hopping to t y = 0.1 in order to effectively be in the 2D 
regime. But since the ratio AE/t y remains close to its 

V = value, the accuracy should not change. That is, 
we expect the error in the interacting Green's function 
Sg ~ 10~ 4 and the error in the interacting density-density 
correlation function SC « 10~ 5 . This gives us a high de- 
gree of confidence in analyzing the properties of the 2D 
interacting system. 

In the presence of t y , the departure from ID behavior 
which is characterized by strong oscillations of n(x) can 
be seen in Fig[T] In the ID systems these oscillations are 
present even in the bulk. For the 2D system, they vanish 
in the bulk, the density becomes uniform as in the case 

V = 0. This is an indication that the dramatic departure 
from the free electron gas seen in the interacting ID chain 
is strongly reduced by t y . This is consistent with the 
relevance of t y or the irrelevance of V (in 2D) found in 
perturbative RG. However, a crossover from a LL to a FL 
would be less apparent in g(x) as seen in FigJ^c). This is 
because in ID, a varies only from when V = to 0.25 
at the QCP V = 2t x . At the same time the exponent 
K = tv /2arccos(—V/2t x ), varies from 1 to 0.5. It is thus 
more favorable to use C2d{x) to analyze the crossover. 



FIG. 4: Ground state energy E gs with respect to the energy 
of disconnected chains for V ■ 
system. 



C2d{x), shown in FiglJJd), has a faster decay than in ID. 
Unfortunately, the actual asymptotic behavior of Cid ix) 
is masked by the presence of the remnant of the Friedel 
oscillations at the ends. Unlike the ID situation, in 2D 
they are 7r-dephased with Cid[x). Hence in order to 
access the asymptotic behavior in 2D, we reverted to even 
L x with odd N ex , for which the Friedel oscillations are 
found to be less severe as shown in Figj3] 

Once in the 2D regime, in order to extract reliable 
correlation exponents, it is important to see how finite 
size effects affect the decay of correlation functions. In 
FigEKa), we show C 2D {x) for 66 x 67, 102 x 103, and 
130 x 131 systems for t y = 0.054. 5E in these systems is 
respectively 0.1273, 0.0675, and 0.0531. It can be drawn 
from the behavior of the 102 x 103 and 130 x 131 systems 
that in the regime t y > SE, aside from edges effects, finite 
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FIG. 5: Linear fit of C(x) for a 130 x 131 systems: in all cases 
ID DMRG (circles), exact (crosses), 2D DMRG (squares), 
exact (pluses); (a) V = 0, t v = 0.025, (b) V = 0.75, t y = 
0.038; (c) V = 1.5, t y = 0.054; _(d) V = 2, t y = 0.073. These 
values of t y are chosen so that t y m 1. 
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FIG. 6: K as function of V: ID DMRG (circles), ID exact 
(crosses), 2D DMRG (squares). 



size effects do not significantly affect the decay of the 
correlation functions in the 2D regime. We thus believe 
that the exponents of C2D (%) that we obtained are very 
close to their value in the thermodynamic limit. Even if 
edge effects are less dramatic when an odd N ex is chosen 
in an even L x system, they nevertheless strongly affect 
the extraction of the correlation exponent. In Figj3jb) 
we show the range of the data used for the extraction of 
K. We arbitrary set x > 10 from the origin and from the 
edge. For this choice it can be seen in FigJ3fb) that we 
are far enough from the upturn of C2d(x) caused by the 
edge. 

Our calculations are made for small values of t y , as 
pointed out earlier, it is expected that when t y « 5E, 
the 2D regime is reached. In order to make this argu- 
ment more precise, we computed the ground-state energy 
E gs = E{ty) — E(t y = 0) as function of t y . Two regimes, 
shown in FigQ] are observed. When t y < pSE, E gs re- 
mains nearly equal to the energy of disconnected chains. 
It would be expected that in this regime, the system will 
essentially have a ID behavior. But when t y > p5E, the 
system gains energy with respect to disconnected chains. 
This is an indication that the system has entered the 2D 
regime. The typical value of p is about 0.5 for the values 
of V we investigated. 

Since we wish to analyze the effects of V on the 2D 
system, we must be careful to avoid a spurious 2D to ID 
crossover which is related to finite size effects. This oc- 
curs when we increase V. Starting at a relatively small 
value of t y and V for which the condition SE < t y is sat- 
isfied for a given size, if we increase V and keep t y con- 
stant, we can reach a regime where t y <C SE. Hence, we 
artificially enter an ID regime. This is clearly a spurious 
effect due to the finite size of the system. In the ther- 
modynamic limit, 5E — > 0, hence, this situation never 
occurs for any finite t y . This problem can be avoided 
by fixing the ratio t y = Jir instead of t y . This means 
that we compensate for the variation of 5E induced by 

V by increasing t y so that we keep the system in the 2D 
regime. 

In Fig[5]we show the decay of Cid(x) and C 2 d(x) for 

V = 0, 0.75, 1.5, 2 in a 130 x 131 system. t y is adjusted 
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FIG. 7: C(x) = x 2 log(x)C(x) relative to C(x = 11) as func- 
tion of x~ 2 for V = 0,0.25,0.5,0.75,1,1.25,1.5,1.75,2 
(from bottom to top); t y ~ 1 in all cases; t y — 
0.025, 0.03, 0.035, 0.038, 0.042, 0.048, 0.054, 0.067, 0.073 re- 
spectively from bottom to top. 



so that t v ~ 1. Therefore, we are in the 2D regime of 
the model. It is important to stress that, if in that case 
a ID like behavior is observed, this would be a genuine 
thermodynamic behavior of the system induced by the 
interactions. FigJSJa) shows that the V = case in ID 
and 2D are consistent with 1/x 2 decay. In FigJSJb), it can 
be seen that the ID data deviate from the 1/x 2 decay. 
On the other hand, the 2D behavior remains similar to 
the V = case. This is consistent with the LL nature of 
the ID system and predictions of a FL ground state in 
2D for mild interactions. When V is further increased we 
find in Figj5jc),(d) that the 2D results deviate from the 
FL behavior as well. Since there is no evidence of CDW 
LRO, this suggests the existence of an unconventional 
metallic state in 2D. A non-linear fit to these data yielded 
K which is displayed in Fig|6] Non-linear fit are known to 
yield many different solutions depending on the starting 
point. To avoid this problem, we first computed the ID 
exponents by fixing the search range in the interval [0, 1]. 
The computed DMRG exponent shown in Fig.© was 
generally in very good agreement with the exact result. 
Surprisingly, the larger discrepancy was observed at V — 
0. These ID exponents were later used as the input in 
the 2D search. The result shows a strong renormalization 
of K from its LL value towards its FL value for V % t x . 
Then, it enters a non-FL regime with K < 1 when t x < 
V < 2t x . 

In Fig. ((TJ) in order to avoid the uncertainty related to 
the fit, the data on C(x) were directly analyzed by study- 
ing C(x) = x 2 log(x)C(x), we added a factor log(x) to the 
FL x~ 2 decay to avoid a maximum that occurs in x 2 C{x) 
at large x and V > 1.25. We believe that this maximum 
is due to logarithmic corrections in the vicinity of the 
QCP. The factor log(x) did not qualitatively modify the 

behavior of C(x) when V < 1. The plot of -£^r for val- 

C(ll) 

ues of t v « 1, shows an evolution from V = to V = 2. 
In the small V regime, C(x) is nearly flat. This is consis- 
tent with the FL physics. There is a downturn at large 
x which is probably due to the influence of the edges. 
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FIG. 8: Spectral Weight function A{Uf , oj) for 64x64 systems 
respectively for increasing height for V = 0.75: t y = 0.005 , 
t y = 0.02, ty = 0.03, and t y = 0.0395 (dotted lines), t y = 0.04, 
t y = 0.05, and t y = 0.06 (full lines). 

In the vicinity of the QCP C(x) increases steadyly with 
increasing x. This clearly proves that K is smaller than 
its FL value in this regime. It should be noted that the 
smallnest of t y implies that the 2D QCP remains very 
close to its ID counterpart. Thus even if V = 2 is not 
exactly at the 2D QCP, it lies very close to it. C(x) 
evolves between these two limits as V increases. This 
shows that the strong V regime clearly departs from the 
FL picture. 



V. LOW FREQUENCY SPECTRUM 

A more direct information on the presence or lack 
of quasi-particles is given by the spectral weight func- 
tion, A(k,w) = E n (*nl4l*o)5(w - Cno), where * is 
the ground-state wave-function, \t„ are the excited state 
wave functions with N e + 1 electrons, and £„o are the ex- 
citation energies between the levels and n. The scope 
of finite frequency study will necessarily be limited to 
very low frequencies. The multiple RG steps in ID and 
2D have truncated out most of the Hilbert space of the 
system. We are left with a very tiny fraction of the total 
number of eigenfunctions and eigenvalues. The essential 
goal of this section is to show that A(k, u) near oj = is 
consistent with our conclusions on C%d{x). 

The DMRG can yield the low oj behavior of A(k, oj) 
by targeting lowest states of sectors with N e and N e + 1 
electrons. The low energy spectrum is then obtained by 
diagonalizing the reduced superblock of size ms2 x ms2 
made of the two external blocks. As for the density- 
density correlation, we restricted ourselves to the central 
chain and used the following approximation, A[k x ,oS) = 

E„(*nl4j*o)e"^^/^F, with e = 10- 2 . This 
means that we are not exactly at the Fermi point of the 
2D systems. In presence of t y , the Fermi point along 
k y = is (kp +Skp, 0), with Skp = Arccos(—ty /tx) — kp . 
Since t y /t x <C 1, Skp is very small. The ID Fermi point 
remains very close to the 2D Fermi surface. 

As in the case of C2d(x), ift y <C pSE, A(k, oj), the sys- 
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FIG. 9: Spectral Weight function AQzf, oj) for 64x64 systems: 
(a) V = 0, t y = 0.05; (b) V = 0.75, t y = 0.05; (c) V = 1.5, 
ty = 0.075; (d) V = 2,t v = 0.15. 



tern will not show the 2D behavior. This is seen in spectra 
displayed in Fig(8]for V = 0.75 and L x x L y — 64 x 65. 
When t y < 0.04 (this is the point where a cusp is seen in 
E gs ), there is a pseudogap in A(k, lu). This pseudogap is 
a ID finite size effect. And as soon as the t y > 0.04, i.e 
the system enters the 2D regime following the criterion on 
E gs , a quasi-particle peak appears in A(k, lo). Hence, to 
some extent if L is large enough, increasing t y is equiva- 
lent to increasing the size of the system. This transition 
is very sharp. This is consistent the cusp observed in 
E gs . In the 2D regimes the peak becomes sharper when 
t y is increased. This was observed for 0.04 < t y < 0.075. 
When t y > 0.075, the peak was less sharp. This is due to 
the fact that for this value of t y , the condition t y <C AE 
was no longer satisfied, hence we started loosing accu- 
racy. 

The spectra shown in Fig [5] are consistent with the 
conclusions drawn from C(x). When V < t x , there is 
a well defined quasi-particle peak at the Fermi energy for 
k x = kp. But when t x ^ V ^ 2t Xl the peak shifts away 
from the Fermi energy; A(kp,tu) displays a quantum- 
fluctuation induced pseudogap which is a precursor of 
the CDW gap. 

The pseudogap observed between 1 < V < 2 exists 
for a finite range of V. This is for instance illustrated 
in Fig. (fit)]) where A(kp,w) is shown for t y = 0.075 and 
t y = 0.1. A pseudogap exists for both values of t y . For 
ty = 0.1, there is a shift of the spectral weigth towards 
lower oj. This is consistent with the possible crossover 
towards a FL at higher t y . 
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FIG. 10: Spectral Weight function A(k F ,uj) for 64 x 64 sys- 
tems for V = 1.5: t y = 0.075 (full line), t y = 0.1 (dotted 
line). 

VI. DISCUSSION AND CONCLUSION 

In this work, we have shown that the two-step DMRG 
is a very useful tool for the study of quasi-lD models. 
This was illustrated in a spinless fermion model. We 
were able to obtain a very good accuracy on systems as 
large as L x [L + 1) = 130 x 131. This is out of the reach 
of the conventional DMRG algorithm 21 . The essential 
difference with the coventional DMRG is the separation 
of the two energy scales of the Hamiltonian. The diffi- 
culty that arised in the interpretation of data was due not 
to accuracy but rather to the effects of open boundary 
conditions. In principle this could be avoided if periodic 
boundary conditions are used. However, we find that 
even periodic boundary conditions were not free of prob- 
lem. When they are applied, the ground state is two-fold 
degenerate for even N e . This means that in constructing 



the ID Hamiltonian, additional effort should be made to 
keep tract of these multiplets. The computed low energy 
Hamiltonian for a single chain is then less accurate than 
with open boundary conditions. 

The analysis of the correlation functions and of the 
low-energy spectral weight function revealed that a FL 
behavior occurs for small interactions. However, when 
the values of the interaction are close to the ID quantum 
critical point, V = 2t x , the system behavior departs from 
that of a FL. A conservative view would be to infer this 
discrepancy to finite size effects. We did our best to 
disprove this possibility by showing that the decay of 
C(x) is essentially identical for 102 x 103 and 130 x 131 
systems for V = 1.5. Given the nearly size independence 
of K between 102 x 103 and 130 x 131, it is unlikely 
that the value K = 0.83 at V = 1.5 and t w 1 would 
significantly change to reach the FL value K = 1 in the 
thermodynamic limit. 

These results could be explained by the unconventional 
ideas raised by Anderson and coworker s 17 ! 18 . When < 
V ^ 1, the LL would be unstable against any small t y . 
The resulting state is a FL. However, when 1 < V < 2, 
though t y is a relevant perturbation, there are no quasi- 
particles until t y exceeds a certain treshold. 
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